Phase Fluctuations and Spectral Properties of Underdoped Cuprates 
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We consider the effect of ciassicai phase fluctuations on the quasiparticie spectra of underdoped 
high-Tc cuprate superconductors in the pseudogap regime above T c . We show that photoemission 
and tunneling spectroscopy data are well accounted for by a simple model in which mean field d-wave 
quasiparticles are semiclasically coupled to supercurrents induced by fluctuating unbound vortex- 
antivortex pairs. We argue that the data imply that transverse phase fluctuations are important at 
temperatures above T c , while longitudinal fluctuations are unimportant at all temperatures. 
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I. INTRODUCTION. 

It is now a well established experimental fact that the 
underdoped cuprate superconductors exhibit a "pseudo- 
gap" behavior above the superconducting critical tem- 
perature T c , which is characterized by a vanishing su- 
pcrfluid density as obtained by transport measurements, 
but persistence of a gap in the quasiparticie excitation 
spectrum as measured by various spectroscopies H . Re- 
cent angle resolved photoemission (ARPES) and 
scanning tunneling spectroscopy (STS) |(^,|7j experiments 
indicate that in the underdoped cuprates the gap evolves 
smoothly as the temperature is increased through T c . Al- 
though the gap begins to fill in and the sharp quasiparti- 
cie peaks are lost above T c , the position of the gap edge 
changes only slightly. This is in sharp contrast to the 
behavior of overdoped cuprates and conventional super- 
conductors where the gap closes at T = T c . It is further 
found that the gap increases as the doping concentration 
is reduced from its optimum value, while at the same 
time T c decreases. This results in highly anomalous ratios 
2A/fcsT c which were reported to attain values of 12 or 
more in Eh^S^CaC^Os+a (BiSCCO), compared to the 
weak coupling BCS value of 3.54. ARPES results also 
indicate that the angular dependence of the gap func- 
tion on the Fermi surface, which below T c follows the 
Ak = Adcos(20) shape expected for a d x 2_ y 2 order pa- 
rameter, develops extended gapless regions around the 
nodes above T c whose size increases with T. Although a 
number of important experimental issues remains to be 
settled, such as the temperature at which the pseudogap 
closes and its possible persistence in optimally and even 
overdoped cuprates suggested by the recent STS results 
[^J, the basic picture of a superconducting quasiparticie 
gap persisting over a wide range of T > T c in underdoped 
materials is well established. 

Many theoretical concepts, including spin fluctuations 
[ |L condensation of preformed pairs ||, SO(5) symmetry 
[ 10| and spin-charge separation pj| , have been invoked to 
explain the pseudogap behavior. In this paper we study 
the implications of a scenario put forward by Emery and 
Kivelson jl2p w ho, following earlier work of Uemura and 
co-workers proposed that the underdoped material 
above T c is in a state with a non-zero local amplitude of 



superconducting pairing, but is not truly superconduct- 
ing due to thermal fluctuations in the phase of the order 
parameter. Within such a scenario the transition at T c 
is of the Kosterlitz-Thouless (KT) type, slightly rounded 
by the weak coupling between the copper-oxygen planes 
along the c-axis. In a strictly 2D system the KT transi- 
tion is associated with proliferation of unbound vortex- 
antivortex pairs. Weak coupling between the planes leads 
to correlated motion between vortices in adjacent planes 
which form 3D vortex loops close to the critical temper- 
ature. The transition to the disordered phase is then 
characterized by the appearance of vortex loops with ar- 
bitrarily large radii [E3. 

In the present paper we study the spectral properties of 
a superconductor in the incoherent state above the phase 
disordering transition but below the mean field transition 
at which the local gap forms. We find that fluctuating 
currents arising from unbound vortex-antivortex pairs 
can contribute significantly to the ARPES and STS line- 
shape broadening in the low-energy region of the spec- 
trum. By analyzing the experimental data we estimate 
the strength of these superconducting phase fluctuations 
and we deduce the vortex core energy. 

We model the underdoped cuprate superconductor as 
a set of independent 2D superconducting layers, each un- 
dergoing a KT transition at a temperature Tkt which we 
identify with the superconducting critical temperature 
T c . The weak interplane coupling, which we neglect, will 
affect the very long lengthscale physics (changing e.g. 
the universality class of the transition from KT to 3D 
XY), but should not affect the shorter lengthscale fluc- 
tuations which contribute to the electron spectral func- 
tions of interest here. The disordered state above Tkt 
can be thought of as a "soup" of fluctuating vortices with 
positive and negative topological charges and with total 
vorticity constrained to zero. Each of these vortices is 
surrounded by a circulating supercurrent which decays 
as 1/r with the distance from the core. Such super- 
currents, within a semiclassical approximation, lead to 
a Dopplcr-shifted local quasiparticie excitation spectrum 
of the form Ofl 



E k = E» + hk-v s (r), (1) 
where v s (r) is the local superfluid velocity and E^ — 
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\J e\ + |Ak| 2 is the usual BCS spectrum. The change 
in the local excitation spectrum will affect the spectral 
properties of the superconductor in that the physically 
relevant spectral function must be averaged over the po- 
sitions of fluctuating vortices. 

This effect will be particularly pronounced in a d-wave 
superconductor since Eq. (Q) implies formation of a re- 
gion on the Fermi surface with E^ < around a nodal 
point for arbitrarily small v s (r). Physically this corre- 
sponds to a region of gapless excitations on the Fermi 
surface which leads to a finite density of states (DOS) at 
the Fermi level. As first discussed by Volovik |17|] , a sim- 
ilar situation arises in the mixed state of a ei-wave super- 
conductor where the superflow around the field-induced 
vortices leads to the residual DOS proportional to \/~H. 
This unusual field dependence arises because the dis- 
tance between vortices in the vortex lattice d v ~ H^ 1 ^ 2 
and the average superfluid velocity projected onto a gap 
node direction is proportional to d~ . At low T and high 
field this implies a ~ T\[H contribution to the electronic 
specific heat which was indeed observed in the measure- 
ments on YBa 2 Cu30g.95 sing le crystals In the 
present case, instead of a regular Abrikosov lattice of 
field-induced vortices, we consider a fluctuating plasma 
of thermally induced vortices and antivortices. The es- 
sential physics however remains the same. 

II. THEORY: QUASIPARTICLE EXCITATIONS 
COUPLED TO PHASE FLUCTUATIONS. 

We shall be interested in how the supercurrents in- 
duced by phase fluctuations affect the spectral function 
of a superconductor, A(k, uS) = — 7r _1 Im(?(k, uS), which 
may be measured by ARPES and STS experiments. Here 
0(k, ui) is the diagonal part of the full superconduct- 
ing Green's function which solves the Gorkov equations 
for a d-wave superconductor, given in the Appendix. In 
the mean field approximation (neglecting, among other 
things, phase fluctuations) the diagonal Green's function 
may be written as 

A 2 

Go 1 (k,uj)=uj-e k + iT 1 Js-, (2) 

W + £k 

where, following ||, we have added to the usual mean- 
field solution a single particle scattering rate Note 
that this form of the scattering rate in Qq constitutes 
an non-trivial assumption. It is not pairbreaking, in the 
sense that it is ineffective at small uj, e^; i.e. in the 
region ui < E^ ~ \Ak + A 2 .. By contrast in a d-wave 
superconductor, a conventional scattering rate enters via 
the replacement to — * ui + iT, leading to a broadening 
which is effective even at low uj, e^. As shown by Norman 
et al. H the form given in Eq. (|J) agrees with the ARPES 
data at T < T c . We demonstrate below that it also agrees 
with STS. 



At T > T c Norman et al. f| showed that addi- 
tional pairbreaking scattering is needed to account for the 
ARPES data, which they modeled phenomenologically 
by introducing another scattering rate Tq ^ Ti, mak- 
ing a replacement ui — ► ui + iT$ in the last term of Eq. 
(g). They suggested that r could arise from exchange 
of pair fluctuations; we find, by explicitly evaluating the 
corresponding propagator [p0| , that this proposed mech- 
anism does not account for the observed magnitude of 
Tq. This conclusion is supported by the results of Vilk 
and Tremblay Q. 

We now discuss what we believe to be a more likely 
source of the pairbreaking scattering, namely supercur- 
rents induced by phase fluctuations. In order to deter- 
mine how Go is changed in the presence of superflow it 
is useful to recall the origin of the energy shift in Eq.(|J). 
This can be derived |l^ , |l6f| by assuming a state of uni- 
form superflow with v s = 7iq/m induced by an order 
parameter of the form e 2iq r Ak. By solving the appropri- 
ate set of Bogoliubov-de Gennes equations and retaining 
only terms to linear order in q, one finds that the en- 
ergy is modified as indicated in (jl]) while the coherence 
factors are to the same order unchanged. This result is 
then semiclassically extended to non-uniform situations 
by assuming slow spatial variations of v s (r). 

One can follow this exact procedure and solve the ap- 
propriate Gorkov equations for Qq in the presence of su- 
perflow. One finds (see Appendix) the following intu- 
itively plausible result which is exact for uniform flow up 
to terms linear in q : 

£ q (k,w) =</ (k-q,w-i7), (3) 

where rj = hvpQc) ■ q ~ Kk ■ v s . Here vp(k) = 
(dek/d\i)k=k F — Kkp/m is the Fermi velocity and the 
last equality holds when the Fermi surface is approxi- 
mately isotropic. In the following we shall assume that 
Eq. (||) can be applied locally when v s (r) varies slowly in 
space. Applying the above prescription to (^) one finds, 
again to the leading order in q, 

g~ l (k, W ) = U £ k + iT x ( ^ k ~° 9 2 , (4) 

where £ = VA(k) • q with VA(k) = (dAk/dk)k=k F ■ One 
can easily estimate v^/vp ~ (fo^f) 1 ~ A^/e^? which is 
typically a small number in superconductor. We there- 
fore expect that Q <C rj. A more detailed numerical analy- 
sis indeed shows that, as long as A^/e^ is small compared 
to unity, the effect of £ on the spectral lineshape is neg- 
ligible compared to that of r/, and will be dropped in the 
following. 

A typical experimentally measured quantity, such as 
the ARPES or STS lineshape, will provide information 
on Qq averaged over the phase fluctuations. Thus, we 
need to evaluate 

£ q (k,w)= / d V P( V )gq(k,u;), (5) 
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where P is the probability distribution of rj given by 

Pfa) = <<*[»? -ftk-v.(r)]). (6) 

The angular brackets indicate thermodynamic averaging 
over the phase fluctuations in the ensemble specified by 
the 2D XY Hamiltonian El 



T~Lxy 



Hx)' 



\{r)d 2 r, 



(7) 



where v s (r) is understood to contain both longitudinal 
(spin wave like) and transverse (vortex like) excitations. 
V = Vo/hsT is a dimensionless coupling constant and Vq 
is related to the superfluid density n s by Vq = fi 2 n s /Am. 
In the nearest neighbor XY model, UbTkt — 0.9Vo [ ji2| ; 
more generally k B T KT = (tt/2)V q (T = T KT ) |§. Lon- 
gitudinal phase fluctuations result in the spatial modu- 
lation of charge density and will be therefore suppressed 
by Coulomb interaction at long wavelengths. This inter- 
action is not explicitly included in the XY Hamiltonian 
(fjj) but we return to it shortly. 

The last term in Eq. (^) can be thought of as a super- 
conducting self energy E s (w,k). Eqs. (||,|]) then imply 
that the primary effect of the phase fluctuations is to 
smear the functional dependence of S s (w,k) on the en- 
ergy variable, broadening the spectral lineshape. A more 
detailed analysis shows that 77 acts primarily to fill in the 
gap, in a way similar to the inverse pair lifetime Tq in- 
troduced by phenomenological considerations in Ref. || . 
r 1; on the other hand, does not affect the lineshape at 
low energies: notice that (?o(kp,cj = 0) = for any Fi. 

We now give a quantitative description of this broad- 
ening by explicitly evaluating P(r)) and the resulting line- 
shapes as a function of temperature. Making use of the 
identity S(x) = (27r) _1 J dse lsx Eq.(^) becomes 



(8) 



To the leading order in cumulant expansion one can write 
( e -»lik.v.0r)} = exp {-i s 2h 2 k a k p (vyP)}, where summa- 
tion over repeated indices is implied and the spatial vari- 
able has been suppressed. This statement becomes exact 
when the transverse fluctuations can be represented by 
Gaussian degrees of freedom, as is done in the Debye- 
Hiickel approximation employed below. The s-integral 
in (@) can now be explicitly carried out, yielding a Gaus- 
sian distribution 



P(ri) = (2ttW) 



-i/2 e -i7w 



(9) 



with W = h 2 k a kf3(vfvP) . We have thus reduced the 
problem of finding the probability distribution P to eval- 
uation of a correlator (vfv^). For a 2D system described 
by the Hamiltonian (7p, this correlator has been consid- 
ered by Halperin and Nelson p7| . They found, using a 
Debye-Hiickel approximation valid for for T well above 
Tkt, 



(«?(P)«?(-P)> 



2m 2 V 



PaP/3 . $a(3 ~ PaPp/p 2 



l + c 3 £ 2 p 2 /4ir 2 V 



(10) 

where v s (p) = J dre~ lp r v s (r). The first term in brack- 
ets comes from the longitudinal and the second from the 
transverse fluctuations. C3 = (2E C — 7r 2 Vb) /ksT is a di- 
mensionless quantity related to the density of vortices, 
E c is the vortex core energy and £ c is the core cutoff. 

Upon averaging over fluctuations the translational in- 
variance is restored and we may evaluate the real space 
correlator at r = 0: 



«(0Vf (0)) = 



1 



{2-Kf 



A«(p)«f(-p)>. (11) 



Explicit integration finally yields, for |k| = fcp, 

V 



w 



8V 



1 



In 1 



C3 



£3 
V 



(12) 



The short wavelength divergence in (|llj) has been cut 
off at q c = 2tt/^ c and we used the BCS relation £ c = 
Hvf /7rArf. The first term in brackets comes from lon- 
gitudinal fluctuations and would be suppressed in re- 
alistic models in which the Coulomb interaction is im- 
portant. The second term, a v = (V/c3)ln(l -I- C3/V), 
comes from transverse fluctuations due to vortices, and 
is always positive and smaller then 1 (this follows since 
C3 > is required by the stability of the system). To the 
extent that E c is independent of temperature the ratio 
C3/V, and therefore a v , is T-independent. The primary 
T-dependence of W therefore comes from V in the pref- 
actor. After expressing Vq in terms of Tkt, W may be 
written as 

W~ 3.48(1 + a v )(T/T KT )A 2 . (13) 

Eqs. (||) and (^|) describe the effect of classical 
phase fluctuations on the spectral function of a supercon- 
ductor. From the knowledge of such a spectral function 
one can compute the respective ARPES and STS line- 
shapes, extract the parameter W, and compare it with 
the prediction given by Eq. ( |l3| ) . This will be done in the 
next section, but we first discuss the validity and some 
qualitative aspects of the results presented above. 

Our results depend on three assumptions; that the 
electron Greens function may be calculated using semi- 
classical methods, that the phase fluctuations are quasi- 
static, and that there is a sufficiently wide KT temper- 
ature regime in which reasonably well defined vortex- 
antivortex plasma exists. The first of these assumptions 
applies when the coherence length is sufficiently larger 
than the inter electron spacing (i.e. kpt; 3> 1) and is the 
same assumption as underlies Volovik's prediction jl7] of 
a ~ T\fH dependence of the specific heat in the mixed 
state of a d-wave superconductor. As this behavior is ob- 
served in YB2CU3O7 jl8|,[l9|], we believe that this assump- 
tion is well justified. The second assumption, of quasi- 
static phase fluctuations, has two parts. The transverse 
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fluctuations come from vortices, so for them the essential 
assumption is that vortices move slowly compared to elec- 
trons. This is justified by Bardeen-Stephen results, which 
imply that vortex motion is overdamped and thus diffu- 
sive, and so surely slower than the ballistic motion of the 
electrons. For longitudinal fluctuations the situation is 
less clear. If the hypothesis of Emery and Kivelson, that 
they are classical (i.e. that at q ~ we have ui q <C T) 
fL2[ is accepted then Eq. (|l^) applies. However, as we 
shall see, the data contradict this. A more likely scenario 
is that Coulomb interaction pushes the longitudinal fluc- 
tuations up to the plasma frequency, in which case the 
coupling to electrons is very weak and one should simply 
remove the longitudinal fluctuations from the theory. We 
shall see that the data are consistent with this picture. If 
(for some as yet unknown reason) the longitudinal fluc- 
tuations are collective modes with a velocity of the order 
of the Fermi velocity, then our results do not apply. The 
final assumption, of a wide temperature regime between 
the mean- field and KT transitions, is the most difficult to 
justify, except on empirical grounds. This hypothesis was 
proposed in Ref. |L2j and our results are consistent with 
it. The theoretical justification for a clean system must 
involve proximity to a Mott insulating phase, which sup- 
presses the superfluid stiffness and hence Tkt, but does 
not suppress the pairing. A detailed theoretical treat- 
ment in two dimensions has not been given. 

As mentioned above, W given by Eq. ( |l3|) describes 
both longitudinal and transverse fluctuations of the phase 
and can be written accordingly as W = Wl + Wt, where 



W L a 3.48(r/r KT )Ad, 
W T =s 3A8a v (T/T KT )A 2 d , 



(14) 
(15) 



The expression for Wt is valid at temperatures well above 
Tkt, where all pairs can be thought of as unbound and 
thermal energy dominates over the intervortex interac- 
tion p^j . Well below Tkt, on the other hand, we ex- 
pect Wt = since the vortices appear only in tightly 
bound pairs which contribute negligible supercurrent be- 
yond the lengthscale set by the pair size and the pair 
density is exponentially small ~ e - 2E ^/ k sT _ gy numer- 
ically integrating the appropriate scaling relations pjj ] 
one could in fact obtain Wt at all temperatures. How- 
ever such level of detail is beyond the scope of this paper 
and we shall confine ourselves to the limiting cases stated 
above and note that Wt has nonsingular monotonic be- 
havior across Tkt- 

The expression (|l4|) for Wl is expected to be valid 
down to low temperatures, provided quantum fluctua- 
tions and the Coulomb interaction can be neglected. In 
a (i-wave superconductor the temperature dependence of 
Wl will be modified by the T-linear temperature depen- 
dence of the superfluid density [28) n s ~ A _2 (T) which 
enters the definition of Vo in the XY Hamiltonian (J?]). 
It is interesting to note that at temperatures below T C1 
say at T = T c /2, Eq. @ implies Wl ^ A|, i.e. large 
broadening of the spectral function by longitudinal fluc- 



tuations. Such a large broadening, comparable to the 
gap itself, would completely obliterate any signature of 
the gap in the excitation spectrum. Clearly, this is not 
observed experimentally ^]|§|6|. As shown below and in 
Ref. [||, experimental data are consistent with W = 
below T c . We must therefore conclude that longitudinal 
fluctuations are strongly suppressed by the Coulomb in- 
teraction as suggested in Ref. . On the same grounds 
we may argue that the observed linear temperature de- 
pendence of the magnetic penetration depth |2£| is not 
due to phase fluctuations, as suggested by some authors 
fi0| , [3l| , but due to thermally excited quasiparticles in the 
nodes of the d-wave gap Q. Transverse fluctuations, on 
the other hand, result in no net charge displacement and 
are thus unaffected by Coulomb interaction. 



III. COMPARISON TO EXPERIMENTAL DATA 

A. Photoemission 

As shown by Norman and co-workers || , a particularly 
simple relationship exists between the spectral function 
of a superconductor and a symmetrized ARPES spectrum 
at the Fermi surface, 



I s (uj) cx A(k F ,u), 



(16) 



where Is(w) — I(u>) + I(—u>) and I{u>) is the measured 
ARPES lineshape. The advantage of this symmetrized 
representation is that the thermal broadening due to the 
Fermi functions is automatically subtracted out. The 
proportionality ( |l6| ) is expected to hold for |w| less then 
few tens of meV and remains valid in the presence of a 
symmetric energy resolution function. We first discuss 
qualitative features of our predicted lineshape and the 
experimental data, and then present results of a detailed 
fitting procedure. 

The symmetrized experimental ARPES lineshapes for 
the 83K underdoped BiSCCO sample of Ref. @ for the 
k = kp vector along the (0,0) — > (0, 70 direction at se- 
lected temperatures are shown in Fig. (j](a). Below T c we 
wish to model these lineshapes by the spectral function 
corresponding to Eq. (||), which for k = kp, e^ F = 
becomes 



Is(u>) oc A (k F , u) = - 



Tr^-A^+Tf 



(17) 



As noted in || , this functional form indeed describes the 
data below T c well, after it is convolved in u) with a Gaus- 
sian of width a — 13.5meV representing the estimated 
experimental resolution. Qualitatively, Ak sets the posi- 
tion of the quasiparticle peaks and Ti (together with 0) 
sets their width. Note that for Ti = cr = the above 
lineshape consists of two (^-functions at w = ±Ai<. It is 
thus clear that fairly large values of Ti in Eq. ( |l7| ) are 
needed in order obtain quasiparticle peaks of the correct 
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width. Even for large Ti the theoretical lineshape I(u) 
tends to zero for \u>\ larger then several Ak, while the 
experimental lineshape saturates to a finite value. Un- 
derstanding this large-w background presents a challenge 
for any theory of ARPES in the cuprates and is beyond 
the scope of this paper. Here we focus on the low energy 
region of the spectrum where we may reasonably expect 
the present simple model to be valid. 

For cr <C Ti one can easily estimate is(0) « 
J (ricr 2 /A k ) and I s (A k ) w Io^f 1 . The peak to valley 
ratio 



_ MA k ) 



ls(0) i> 



K_ 

2^2 



(18) 



is independent of the unknown prefactor Iq and can be 
easily extracted from the raw data. Assuming fixed Ak 
one can thus obtain reliable estimates of I\ without per- 
forming detailed fits. In particular, application of Eq. 
( |l8| ) to the data in Fig. E(a) implies that Ti grows by 
about a factor of 6 between T = 14K and 100K. The 
above analysis also implies that at low ui the lineshape 
depends crucially on the experimental resolution a. For 
instance decreasing er by a factor of 2 the ratio k should 
grow by a factor of 4. Confirming this prediction exper- 
imentally would be a valuable test of the present model. 

Above T c Eq. (|l7]) no longer provides a good fit for 
the data. The reason for this is the persistence (up 
to ~ 200K) of a well defined edge-like feature around 
lu = Ak along with a pronounced increase in the low-w 
density of states. This behavior cannot be modeled by 
further increasing T% since the values needed to fix 1(0) 
would rapidly smear the edge. Inclusion of the phase fluc- 
tuations, i.e. finite W in the averaged spectral function 
(||), rectifies this problem. Analytically it is somewhat 
difficult to discuss the combined effect of W, T\ and a 
on the lineshape. Qualitatively one can show that the 
primary effect of increasing W is to "fill in" the gap. 
This is precisely what is needed to describe the ARPES 
data above T c . 

Fig. |l|(a) shows our fit to the symmetrized ARPES 
lineshapes for the underdoped sample using a numeri- 
cal computation of the full spectral function extracted 
from Eq. (^|). Least square fits in which Ak, Ti and 
W were taken as free parameters, were performed after 
convolving the spectral function _A(kp,u;) with experi- 
mental resolution a — 13.5meV. In the low energy region 
\lo\ < 80meV, where the simple model Green's function 
approach with w-independent scattering rate Ti is ex- 
pected to be valid, the fits are excellent for all temper- 
atures. The extracted parameters are displayed in Fig. 
|l|(b). Both Ak and Ti behave in the way expected for an 
underdoped cuprate: the gap is approximately constant 
across T c while the scattering rate rises sharply below 
T c and saturates at higher temperatures. In the present 
model the large increase in Ti is required to wipe out the 
quasiparticle peaks. W also behaves as anticipated from 
the above considerations. At low temperatures W ~ 
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FIG. 1. a) Symmetrized ARPES lineshapes Is(u) at mo- 
mentum near (n, 0) from Ref. j^| (thin line) and our fits to the 
spectral function (thick line) at selected temperatures. The 
dashed line illustrates fit with W = 0. Region of the fit is 
|w| < 80meV. b) Parameters of the fit as a function of T. 



indicating that below T c the phase fluctuations are neg- 
ligible; vortices appear only in tightly bound pairs and 
longitudinal fluctuations are suppressed by the Coulomb 
interaction. Above T c fluctuations become important 
(pairs unbind) and W approaches the T-linear behav- 
ior consistent with ( |l3| ) . We stress here that a good fit to 
the data requires W > above T c . This is illustrated by 
the dashed line in Fig. |l](a) which represents the fit to the 
100K lineshape with W = and the gap value restricted 
to A k < A k (T = 0) §||. We note, however, that in this 
model much of the observed Iow-cj density arises from the 
large value of Ti in combination with the instrumental 
resolution. It is possible (and is indeed suggested by the 
analysis of the STS data in the next section) that the 
large T\ contribution is an artifact, arising from a com- 
bination of experimental resolution in the photoemission 
experiments and inadequacies of our theoretical model. 
We therefore regard the values of W obtained here as 
underestimates. 

From the slope of W(T), assuming that transverse fluc- 
tuations are dominant, we estimate a v ~ 0.009 imply- 
ing the vortex core energy E c /Vo ~ 360. This value is 
much larger than the usual condensation energy in the 
vortex core E c ~ 2Vb Jr|. Within the Debye-Hiickel 
approximation the vortex density can be estimated as 
Pv = ^- 2 (TT/c 3 )(l-a v ) ~ ^- 2 (T/T KT )(V/c 3 ), for a v « 1 
and fceXxT — Vo. This implies, for the parameters ex- 
tracted from ARPES data, p v ~ 4.3 x 10~ 3 £- 2 (T/T KT ). 
It is remarkable that such a small density of vortices leads 
to significant broadening of the lineshape. We should also 
remark that for such a small density of vortices one may 
question the validity of Debye-Hiickel approximation at 
temperatures in consideration. We emphasize, however, 
that only our interpretation of W and in particular the 
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estimates of E c and p v depend on the validity of this ap- 
proximation. Our analysis of the lineshapes is quite gen- 
eral, since we treat W as a free parameter of the model. 

We have performed similar fits for an overdoped 82K 
BiSCCO sample of Ref. §. Our results are consis- 
tent with Ak vanishing close to T c and, within statis- 
tical noise, W — at all temperatures. This indicates 
that phase fluctuations are unimportant in overdoped 
cuprates and the transition is essentially mean-field-like. 

We now consider ARPES data at the Fermi crossing 
close to the gap node direction (it, ir). These indicate ex- 
tended regions of gapless excitations above T c which grow 
in size with temperature [§||]5)- This is not reproduced 
by the simple model we have considered so far. The rea- 
son is that adding the /iv s • k term to the quasiparticlc 
energy [cf. Eq.(||)] effectively depletes the local spectral 
function for v s parallel to k but enhances it by equal 
amount for opposite v s . Upon averaging over all direc- 
tions of v s (for fixed k) the net effect is to broaden the 
mean- field lineshapes as seen in Fig. 0(a). Phase fluctu- 
ations cause no net depletion of the spectral weight near 
the gap nodes. 

In order to account for the ARPES data the form of 
Ak must change. Since the supercurrent flowing around 
individual vortices is pair-breaking, it is in principle pos- 
sible that it will alter the internal structure of the self- 
consistent gap function in addition to usual suppression 
of the order parameter in the core. A similar scenario 
has been proposed by Haas et al. |35| who considered the 
effect of non-magnetic impurities on a d-wave gap func- 
tion. They found that it was possible to construct a gap 
function such that with increasing disorder nodes indeed 
expanded into finite gapless arcs. Intuitively this effect 
can be understood on the grounds that smaller gap near 
the nodes is more susceptible to the pair breaking. In 
the present case pair breaking is caused by supercurrents 
rather than impurities, but the physics remains the same. 

In order to substantiate this idea we have solved the 
self-consistent gap equation for a <i-wave superconduc- 
tor in the presence of uniform superflow. We considered 
a model gap function of the form Ak = A^ cos(2#) + 
A^ cos(6(9) with the appropriately generalized pairing in- 
teraction |Q. We found that, for a system with A^ < 
in the absence of superflow, a transition to A^ > occurs 
when sufficiently strong supercurrent flows in the direc- 
tion close to the nodal vector k„; i.e., when |v s • k n | > 
bAd/hkp, with b ~ 0.7 a model dependent constant. The 
state with A^ > exhibits extended gapless regions (cf. 
Fig. ^) while that with A^ < only the usual point 
nodes. Extrapolating this behavior to the supercurrent 
flowing around the vortex we argue that a region of ex- 
tended gapless excitations may form in the vicinity of 
the core. Such a region would have the shape of a four 
leaf clover (schematically depicted in the inset of Fig. ||) 
and a spatial extent of several £ c . A truly quantitative 
treatment of this effect is complicated because it involves 
sclf-consistently solving the d-wave vortex problem which 
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FIG. 2. Model gap function A(0) = A d cos(26»)+A^ cos(6#) 
with r = A' d / Ad positive (extended nodes) and negative 
(point nodes). Inset illustrates the proposed real space elec- 
tronic structure of the vortex in a d-wave superconductor. 

is a highly non-trivial task |36|. We note, however, that 
recent STS data on vortices in BiSCCO |37j , show a pe- 
culiar pseudogap behavior near the vortex core, which 
may be indicative of a formation of the gapless regions 
around a vortex proposed above. 

Assuming that this picture is correct, it is clear that 
in the vortex-antivortex plasma above T c , upon averag- 
ing over fluctuations, ARPES will detect a gap function 
strongly suppressed for k close to the nodes. Further- 
more, with increasing temperature the volume fraction of 
gapless regions will grow (since the vortex density grows) 
leading to larger gapless areas on the Fermi surface in 
agreement with experimental observation. 

B. Tunneling 

Tunneling conductance, the quantity measured by 
STS, is related to the spectral function of the supercon- 
ductor by 

/oo 
dwf{u-E)Y,\M^)\ 2 A{^), (19) 
-oo k 

where / is the Fermi function and Mk(w) is the tunnel- 
ing matrix element, usually approximated by a constant. 
Tunneling conductance reflects the spectral function av- 
eraged over the entire Brillouin zone and broadened in 
the energy variable by the Fermi function. Thus, unlike 
in the ARPES lineshape function Ig{E), quasiparticle 
peaks in g{E) will be broadened even in the absence of 
scattering and at T = 0. For measurements performed on 
similar samples one would thus expect ARPES spectra 
to be much sharper than STS spectra, at any tempera- 
ture. Fig. ||(a) displays g(E) as measured by STS on 83K 
underdoped BiSSCO sample of Ref. ||. Surprisingly, we 
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FIG. 3. a) Fit to the STS data of Ref. |J at selected tem- 
peratures, b) Parameters of the fit as a function of T. Above 
T c both Ad and Fi are fixed to their values at 84K. 

observe that at lowest available temperatures STS line 
is in fact sharper than the ARPES line. As discussed in 
the previous section the broadening of the ARPES line 
comes exclusively from the scattering and experimental 
resolution. Therefore, inspection of the raw data sug- 
gests that the scattering rate Tf TS will be much smaller 
than r^RPES This conclusion is indeed borne out by 
the detailed fitting procedure carried out below. This is 
a rather surprising result which we discuss more fully in 
the next section. 

As seen from Eq. (|l^) , modeling of the tunneling con- 
ductance requires knowledge of the band structure away 
from the Fermi surface and is therefore somewhat more 
involved than that of ARPES. Nevertheless, we find that 
assuming a simple free electron dispersion with cylindri- 
cal Fermi surface and Ak = AdCOs(26*), provides a rea- 
sonable fit for the data at low temperatures, provided 
that one compensates for the asymmetric background 
conductance and assumes ^-dependent matrix element 
Mit(w) tx |cos(20)|. The latter assumption is motivated 
by band structure calculations Q which indicate that 
tunneling between layers is dominated by the regions on 
the Fermi surface close to the (0, 7r) points. As shown 
in Fig. ||(a), this ^-dependence allows to simultaneously 
account for the sharp quasiparticle peak and the wide 
gap in the STS line. Assuming Mk(w)=const leads to 
broader peaks and sharper V-shaped gap structure, in- 
consistent with experimental data. We further remark 
that our simple model does not (and is not expected 
to) explain the pronounced dip feature appearing in the 
spectrum at higher energies, whose origin is at present 
unknown. Over the range T < T c Tf TS has similar quali- 
tative behavior as r^ RPES [compare Figs. |l](b) and||(b)], 
but remains about one order of magnitude smaller. We 
have attempted to reconcile the two sets of data by con- 
sidering a more realistic band structure and ^-dependent 



scattering rate in Eq. ( |l9|) . However, even under the most 
favorable conditions r^ iS remains factor of 5-6 smaller 
than r£ RPES . 

Above T c the STS line does not have enough features 
to permit a meaningful three parameter fit. In particular 
the sharp gap edge completely disappears above 100K 
which leads to ambiguity in defining A^ from the data. 
Based on our previous finding (from the ARPES data) 
that both Ad and Ti exhibit only weak T-dependence 
above T c , we fix these two parameters at their respec- 
tive values at 84K (40meV and 34meV) and extract the 
temperature dependence of W. This appears to us as a 
reasonable procedure which we further check by perform- 
ing two-parameter fits with Ad held constant or slowly 
decreasing with temperature as implied by Fig. 0(b) . W 
is shown in the inset to Fig. 0(b). Its qualitative behav- 
ior is similar to that of W AT!PES but the amplitude is 
roughly factor of 10 larger. The difference is caused in 
part by the much smaller value of Ti discussed above. 
Part of the discrepancy between t^arpes an( j rySTS can 
presumably be attributed to differences in material and 
experimental uncertainties, as well as the failure of our fit 
to account for the temperature variation of Ti, which ac- 
cording to Fig. |l](b) grows by another factor of 2 between 
80K and 200K. Nevertheless, after accounting for these 
factors, considerable discrepancy remains in place which 
is not understood at present. We estimate a v ~ 0.097 
which implies the vortex core energy E c ~ 22Vo and vor- 
tex density p v ~ 8.7 x 10~ 2 ^ 2 (T/T KT ). The value of 
E c /Vq is still large compared to the conventional estimate 
of the condensation energy in the core E c ~ 2Vq p6[ , but 
is consistent with large core energy E c ~ 26Vo deduced 
from lower critical field measurements of YBa2Cu306.95 
at T = |Q. We note that E c is a cutoff-dependent 
quantity and therefore the precise numerical value quoted 
here should be accepted with that in mind. It is also 
possible that the large value of the ratio E c /Vq is due 
to an unusually small Vq rather than unusually large 
E c . Indeed, the vortex core energy is typically of the 
order of Fermi energy. Estimate of Vq given in Ref. 
fH, 5meV< V < lOmeV, implies E c ~ 0.1 - 0.2eV, in 
reasonable agreement with many theories of underdoped 
cuprates which suggest a Fermi energy of the order of the 
exchange constant J ~ 0.15eV. 

IV. DISCUSSION 

The qualitative behavior of ARPES and STS line- 
shapes in underdoped BiSCCO clearly establishes the 
existence of a scattering mechanism which becomes op- 
erative at T > T c and which acts primarily to fill in 
the gap at low energies. We have shown that trans- 
verse phase fluctuations associated with proliferation of 
unbound vortex-antivortex pairs in the system provide 
a reasonable explanation for this scattering. Our anal- 
ysis also indicates that longitudinal (spin wave) fluctua- 
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tions are almost completely suppressed, above and below 
T c . It has been proposed [|30|j3 1 ] that in high-T c materi- 
als, longitudinal phase fluctuations governed by the XY 
Hamiltonian (Q) are important in that they significantly 
contribute to the observed temperature dependence of 
the magnetic penetration depth [2^]. We have calculated 
the broadening of the spectral function which would be 
caused by these fluctuations, and found it to be much 
greater than the experimental data would permit. We 
therefore conclude that longitudinal fluctuations are sup- 
pressed, perhaps by the Coulomb interaction as suggested 
in Ref. g|. 

Quantitatively there exists considerable discrepancy 
between the parameters describing the ARPES and STS 
lincshapes, in particular the single particle scattering rate 
Ti and phase fluctuation broadening W. Since the dis- 
crepancy is apparent at low temperatures and in over- 
doped cuprates we are led to believe that the problem 
lies primarily in our lack of detailed understanding of 
the lineshapes rather than the physics of phase fluctua- 
tions above T c . The most disturbing is almost an order 
of magnitude difference between r ARPES anc j pSTS f oun( j 
below T c , which is implied directly by the raw data. In 
view of the fact that both measurements pertain to un- 
derdoped BiSCCO crystals with similar critical temper- 
atures, it appears unreasonable to attribute such a large 
discrepancy to the material differences. We speculate 
that the large scattering rate needed to fit the ARPES 
data is an artifact related to our incomplete understand- 
ing of the photoemission process in the superconductor 
which is theoretically not completely understood even in 
simple metals p5f . Tunneling spectroscopy, on the other 
hand, is a technique well established in superconductors. 
We therefore surmise that parameters obtained from STS 
more directly reflect the underlying physics. Indeed 
pSTS ^ g me y a t 4.2K is comparable to the scattering 
rates deduced from transport measurements |39|]40(| on 
underdoped cuprates, and £f TS — 22Vb, although large 
for a conventional superconductor, is perhaps not un- 
reasonable in cuprates p4| . Consequently, ARPES line- 
shapes appear to reflect significant extrinsic broadening 
of unknown origin. The puzzling aspect of this interpre- 
tation is that the additional physics in the ARPES spec- 
tra enters as a multiplicative rather than additive factor 
to the apparent scattering rate; cf. pARPES ^ gpSTS 
over the entire temperature range below T c , in which Ti 
changes by a factor of 6. It is also possible that in the 
cuprates the c-axis tunneling matrix element Mk(w) in- 
troduced in Eq. (19) is itself anomalous. Improving the 
energy resolution a of ARPES could shed some light on 
this issue. As noted below Eq. (|||) the ARPES line- 
shapes are strongly affected by experimental resolution 
at small u>. If the model Green's function (||) is correct, 
a factor of two improvement in a should lead to consid- 
erable decrease in the measured intensity at u> = but 
almost no change in the width or height of the quasipar- 
ticle peaks at \lu\ = A^. 

Finally we note that sizable transverse phase fluctua- 



tions implied by this work will also affect other properties 
of the underdoped systems, such as the electronic spe- 
cific heat, fluctuation diamagnetism and transport. Vor- 
tices existing above T c should also generate local mag- 
netic fields which are zero on average but have a non 
vanishing variance. If such fields could be detected, e.g. 
by muon spin rotation experiment, this would constitute 
a direct evidence for the phase fluctuation model of the 
pseudogap phase. 
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APPENDIX: GORKOV EQUATIONS IN THE 
PRESENCE OF SUPERFLOW 

Real-space Gorkov equations EJ generalized to 
anisotropic superconductors read 

(ui - H e )g(r 1 ,r 2 ;uj) + AJ"+(ri, r 2 ; u>) = S(r 1 - r 2 ), 
(w + HDF+iri, r 2 ; u) + A*0(r 1( r 2 ; u) = 0. 

(Al) 

Here Ti. e = [iV — (e/c)A] 2 /2m — ep is the single electron 
Hamiltonian and A is the gap operator for spin singlet 
superconductivity defined as 

AjF+(n,r 2 ;a;)= / ' d 2 r> A(r 1; r').F + (r', r 2 ; W ). (A2) 

A(ri,r 2 ) is the gap function which is in general a 
nontrivial function of both electron coordinates in the 
anisotropic superconductor. We are interested in the 
state of uniform superflow induced by the gap function 
of the form 



A(r 1 ,r 2 ) = A (r 1 -r 2 )e^( ri+r2 ) 



(A3) 



and A = 0. The easiest way to solve ( |Al| ) for Q is to 
perform a gauge transformation to the gauge where the 
order parameter is real and independent of the center of 
mass coordinate R = (ri + r 2 )/2: 

A(n,r 2 )-> A(r 1 ,r 2 )e-^+ r *) ) 

A -> A - -q, 

e 

g(r U T 2 ;u) -» £(!■!, r 2 ; W ) e -^ ri - r2 >, 
^+(n,r 2 ; W ) -» T + (r ll r 2 ;u J )e^ (r ^K (A4) 
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I t is easy to verify that under such transformation Eqs. 
(Al) remain invariant |4l|| . In the new gauge Gorkov 
equations are manifestly translationaly invariant, i.e. in- 
depend of R. Fourier transforming in the relative coor- 
dinate r = ri — r 2 leads to algebraic equations for Q and 
T + of the form 



(w - e k _ q )£(k, u) + A k ^+(k, u) = 1 
(w + e k+q )^+(k, u) + A k £(k, u) = 0, 

where e k = k 2 /2m — ep. The solution for £/ is 



£k4 



£k-c 



(A5) 



(A6) 



Expanding e k ± q to leading order in q we obtain 



Sq^k,^) = (w - r?) - e k - 
= ^(kjW - 17), 



(w - 77) + e k 



(A7) 



with 77 = vy?(k) • q ~ k • v s . As a final step we transform 
Q back to the original gauge with A = 0. According to 
( A4) this amounts to simply replacing k — > k — q on the 
right hand side of Eq. (A7). We thus obtain the desired 
expression In deriving this result we have assumed 
for simplicity a free particle form of the single electron 
Hamiltonian TL e . Evidently, the calculation remains valid 
for more complicated Hamiltonians. 
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